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By making use of the generalized Landau level representation (GLLR) for the quasiparticle prop¬ 
agator, we study the effect of screening on the properties of the quantum Hall states with integer 
filling factors in graphene. The analysis is performed in the low-energy Dirac model in the mean- 
field approximation, in which the long-range Coulomb interaction is modified by the one-loop static 
screening effects in the presence of a background magnetic field. By utilizing a rather general ansatz 
for the propagator, in which all dynamical parameters are running functions of the Landau level 
index n, we derive a self-consistent set of the Schwinger-Dyson (gap) equations and solve them nu¬ 
merically. The explicit solutions demonstrate that static screening leads to a substantial suppression 
of the gap parameters in the quantum Hall states with a broken t/(4) flavor symmetry. The tem¬ 
perature dependence of the energy gaps is also studied. The corresponding results mimic well the 
temperature dependence of the activation energies measured in experiment. It is also argued that, 
in principle, the Landau level running of the quasiparticle dynamical parameters could be measured 
via optical studies of the integer quantum Hall states. 

PACS numbers: 73.22.Pr, 71.70.Di, 71.70.-d 


I. INTRODUCTION 

As predicted theoretically more than three decades ago 0 , 0 , the low-energy excitations in planar graphite (also 
known as graphene) are described by (2 -|- l)-dimensional, massless Dirac fermions. Interestingly, the spinor structure 
of the corresponding fermion fields has nothing to do with the usual spin of electrons. It is connected with an effective 
“pseudospin” that has its roots in the hexagonal arrangement of carbon atoms in graphene, which can be viewed as a 
superposition of two inequivalent triangular sublattices. At the same time, the usual spin plays a rather passive role 
analogous to an extra “species” degree of freedom. 

The Dirac nature of excitations was confirmed experimentally by the observation of the anomalous integer quantum 
Hall (QH) effect in Refs. ii- At sufficiently low temperatures, the measured Hall conductivity revealed a set of 
well-resolved QH plateaus at the filling factors iz = ±4(fc -|- 1/2), where k = 0,1,2,... is an integer. One of the 
signature properties of the corresponding QH effect is the anomalous shift of 1/2 in the expression for the filling 
factor. This measurement appears to be in perfect agreement with the theoretical predictions of the low-energy Dirac 
theory 0^ • The overall factor of 4 in the expression for iz is also in agreement with the predicted fourfold (spin and 
sublattice-valley) de gen eracy of the Landau levels in the low-energy theory. For reviews of quantum Hall physics in 
graphene see Ref. 0-[ll|. 

The spin and sublattice-valley degeneracy of the Landau levels in the low-energy theory of graphene can be associated 
with an approximate global “flavor” 1/(4) symmetry. Strictly speaking, the symmetry is not exact. There exist several 
small explicit symmetry breaking effects that lift the degeneracy of the Landau levels [l0 . The most obvious among 
them is the Zeeman effect that breaks the t/(4) symmetry down to U^{2) x Ui{2), where Us{2) with s =t,4. is the 
sublattice-valley symmetry of quasiparticles with a fixed spin. In practice, however, the flavor U (4) symmetry could 
be treated as exact because, for any realistic value of the magnetic field, the Zeeman energy is much smaller than 
the Landau energy scale. Moreover, the Zeeman energy, as well as all other explicit symmetry breaking effects are 
negligible not only compared to the Landau energy scale, but even compared to the nonzero (thermal/interaction) 
widths of the individual Landau levels. This explains why it is hard to resolve experimentally the QH plateaus with 
any filling factors other than iz = ±4(fc + 1/2) that correspond to fully filled fourfold-degenerate Landau levels. 

The subsequent experimental studies revealed that, in sufficiently strong magnetic fields, additional QH plateaus at 
iz = 0, ±1, ±3, ±4 can be resolved [isl - l^ . As is clear, the corresponding QH states have fractional fillings of the low- 
lying Landau levels. The implication of this is that the approximate fourfold degeneracy of the Landau levels is truly 
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lifted. Such new QH states could be explained theoretically if the electron-electron interaction triggers spontaneous 
breaking of the flavor t/(4) symmetry. For example, the symmetry breaking C/(4) —>■ C/|(2) x U^{2) by a dynamically 
enhanced Zeeman splitting could potentially explain the origin of the QH states with v = Q and v = ±4. If realized, 
such a scenario would be nothing else, but a textbook example of the QH ferromagnetism (QHF) [2^-[3l| [see also 
Ref. [13]. 

It should be emphasized, however, that QHF is not the only possibility. There exist a number of symmetry break¬ 
ing mechanisms and residual symmetries consistent with the fillin g fa ctors of the additional QH states. One of such 
alternative scenarios utilizes the idea of magnetic catalysis (MC) The corresponding order parameters are 

excitonic^article-hole) condensates responsible for the generation of the Dirac and/or Haldane masses of quasipar¬ 
ticles [H. I36l - l43l |. At the microscopic level, the excitonic condensates corresponds to a char ge d ensity wave (CDW), 
or a valley polarized CDW. Symmetry arguments, as well as direct effective model studies |4(TI . EH suggest that the 
order parameters associated with the MC and QHF scenarios necessarily coexist. Also, in principle, these two could 
reproduce all integer QH plateaus observed experimentally in strong magnetic fields. With that being said, a num¬ 
ber of different types of order parameters are possible It is also fair to note that the precise nature of the 

observed integer QH states is not always unambiguous from theoretical considerations and not always established 
unambiguously in experiment. 

From the experimental point of view, a lot of effort was devoted to revealing the underlying nature of the strongly 
insulating ly = 0 QH state, associated with half-filling of the lowest Landau level [l3l - [^ . The main advances in 
resolving the nature of the QH states were made by applying a tilted ma gne tic field to high-quality graphene devices, 
that were fabricated on a thin hexagonal boron nitride substrate |23l . |24| or encapsulated between two layers of 
hexagonal boron nitride 1^. In the case of = 0 QH state, in particular, a careful analysis of the conductance and 
the bulk density of states [231 - 123 suggests that the state is not spin polarized when the magnetic field is perpendicular 
to the plane of graphene. This is consistent, for example, with both an antiferromagnet and a charge density wave. 
When an in-plane component of the magnetic field increases, the state gradually transforms into a fully polarized 
ferromagnetic state. In the intermediate regime, a canted antiferromagnetic state is presumably realized |44l |. Such 
an interpretation is supported by the observation of the bulk gap that does not close and the edge states that become 
conducting with increasing of the total magnetic field. Considering, however, how heavily the arguments rely on the 
properties of the edge states, the final conclusions should be still accepted with a caution. 

The high-quality graphene devices in Refs. [13 [13 also reveal a large sequence of integer QH states with n > 1. 
Among these, there are two states associated with quarter filling of the lowest Landau level, i.e., v = ±1. The linear 
growth of the energy gaps in these two states as a function of the total magnetic field points to their spin polarized 
nature and the major role played by the Zeeman energy in aligning the ground state. A similar sensitivity to the 
Zeeman energy is also observed for the states associated with half-filling of higher Landau levels, i.e., v = 4n, but not 
for the states with quarter and three-quarter fillings, i.e., v = 4n± 1 [13- As we will see, most of these features are 
reproduced in the model studied here. 

In the mean-field approximation, the role of the long-range Coulomb interaction and Landau level mixing were 
studied in Ref. [13. By utilizing a rather general combination of the MC and QHF order parameters, the corresponding 
study was able to reproduce all observed QH plateaus (i.e., v = 0, ±1, ±3, ±4) as well as to suggest that QH plateaus 
with any integer filling are possible. The symmetry breaking patterns of the solutions obtained in the model with 
the long-range Coulomb interaction appeared to be similar to those in the model with local four-fermion interaction 
in Ref. [4l[. The qualitative effects due to the long-range force were (i) Landau level mixing and (ii) “running” 
of all dynamical parameters (i.e., the wave-function renormalization parameters, the Dirac and Haldane masses) as 
functions of the Landau level index. 

In this paper, we extend the study of the abnormal integer QH effect in the model with the long-range Coulomb 
interaction |48l | by including the screening effects and the effects of nonzero temperature. By making use of a similar 
mean-field approximation that neglects the fluctuations of order parameters, we will derive the full quasiparticle 
propagators in the GLLR formalism for all qualitatively different series of QH states with integer filling factors. (A 
semirigorous justification of the approximation used will be provided in Sec. Ill CI l As expected, the corresponding 
results contain not only the information about the energy gaps at the Fermi level, but also the complete dispersion 
relations of quasiparticles in all Landau levels. The resulting propagators can be used in transport calculations, 
predictions of various emission and absorption rates, etc. 

This paper is organized as follows. In Sec. m we introduce the model, define the quasiparticle propagator, and 
derive the gap equations. In Sec. IHII we present our numerical analysis of the gap equations and classify the main 
types of solutions. We also compare our results with those in the previous studies. The summary of the main results 
is given in Sec. |TVl 
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II. THE MODEL 


Following the same approach as in Ref. [1^ , in this study we will use the language of the low-energy theory for the 
description of the QH effect in monolayer graphene. The low-energy quasiparticle fields are given by the following four- 
component Dirac spinors dtg = {'4’kAs,'4’kbs,'4’k'Bs-, —'4’k'As)-, which combine the Bloch states on the two different 
sublattices {A, B) of the hexagonal graphene lattice in coordinate space. The components labeled by the valley indices 
K and K' correspond to the Bloch states with the momenta from the vicinity of the corresponding inequivalent Dirac 
points [K or K') in the two-dimensional Brillouin zone. The field also carries an additional index that describes 
the spin state, i.e., s =t,4” 


A. Quasiparticle Hamiltonian 


The free quasiparticle Hamiltonian has the following pseudorelativistic form: 


Ho = VF 


d^r'I'(7V2, -I- q^TTy)'!', 


( 1 ) 


where vf ~ 10®m/s is the Fermi velocity, r = {x^y) is the position vector in the plane of graphene, and dt = jg 
the Dirac-conjugated spinor. Note that the sum over the spin states is implicit in Eq. ©■ The canonical momentum 
TT = —i/iV —eA/c includes the vector potential A that corresponds to the magnetic field component Bj^ orthogonal to 
the plane of graphene. The 4x4 Dirac matrices 7 *^ (with = 0,1, 2) are defined as follows: 7 ^' = (g) (r^, ir^, 

where the two sets of Pauli matrices f and r act in the valley {K, K') and sublattice (A, B) spaces, respectively. They 
satisfy the usual anticommutation relations { 7 ^, 7 “^} = , where = diag(l, — 1 , — 1 ), and z/ = 0 , 1 , 2 . 

In order to be able to describe QH states with arbitrary filling factors, we should also allow for a nonzero electron 
chemical potential The latter is introduced via an additional term, —fJ-eJ in the free Hamiltonian in 

Eq. dl]). The long range Coulomb interaction is included in the low-energy Hamiltonian by adding the following term: 

Hc = ^J d2rd2r'4''f(r)4'(r)C/c(r - r')4'^(r')4'(r'), ( 2 ) 

where 17c(r — r') is the coordinate-space Coulomb potential in the presence of a constant magnetic held. It is 
straightforward to check that the resulting Hamiltonian H = Hq + He — fJ-e J is invariant under the flavor 

C/(4) symmetry that combines the transformations in spin and sublattice-valley spaces. Eor the explicit form of the 
symmetry generators see, for example, Ref. @ . Such a havor symmetry is explicitly (although weakly) broken by the 
Zeeman interaction of the quasiparticle spin with the external magnetic held B. The additional interaction term in 
the Hamiltonian is given by J d^r/rsH'I'V^'I', where B = |B|, ys = eh/{2mc) is the Bohr magneton, and is the 
Pauli matrix acting in the spin space. As is easy to check, the Zeeman term breaks the ?7(4) symmetry down to the 
U^{2) xUi{2) symmetry. 

The complete model Hamiltonian, including the Zeeman interaction and the electron chemical potential /Xe, can be 
conveniently rewritten in the following quasirelativistic form: 

id = y [nF( 7 r • 7 ) - ^e7° + dx -I- He- (3) 

This implies that the inverse bare quasiparticle propagator is dehned by the following expression: 

= [{uj + g,e - fJ-BBa^)-f° - VF{Tr ■ j)] S{r-r'). (4) 

Here and below, it is convenient to use a mixed (w, r)-representation. 

Let us note in passing that a certain degree of disorder is always present in real graphene devices and, in fact, 
plays a critical role in the observation of the quantum Hall plateaus. However, in the model analysis below, which 
concentrates primarily on the delocalized quasiparticles states in the bulk of graphene, it is justifiable to ignore 
disorder. Indeed, when the spectrum of the delocalized states is established, a large number of bulk observables (e.g., 
symmetry properties of the QH states, density of states, energy of transitions lines, activation energies, etc.) will be 
predicted without much ambiguity. The corresponding limitation is not so critical also because of the robustness of 
the QH effect that stems from its topological nature. 
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B. Dynamical symmetry breaking and full quasiparticle propagator 

In order to be able to describe various QH states with integer fillings factors, we use a rather general ansatz for the 
full quasiparticle propagator 


G ^(w;r,r') = —z ( 7 ^ 0 ; — z;_FF'*'( 7 r • 7 ) + S’*' (5(r —r'), 


(5) 


where F~^ and E’*' are operator-valued wave-function renormalization and self-energy functions, respectively. For 
simplicity, here we will assume that the propagator is a diagonal 2 x 2 matrix in the spinor space, or in other words 
that the propagator for each of the two spin states looks like that in Eq. As is clear, such a choice of the ansatz 
for the propagator does not allow spin mixing. As a result, one cannot describe certain types of states (e.g., a canted 
antiferromagnetism Hi). Nevertheless, we emphasize that the ansatz in Eq. ([S]) is extremely flexible and can describe 
a large number of states (with very different symmetry properties) for each integer filling factor. In the case of filling 
factor zz = 0 , for example, these include charge and spin density waves, as well as ferromagnetic and antiferromagnetic 
states. These are basically all main options proposed for a configuration with a magnetic field perpendicular to the 
plane of graphene. 

By construction, both F'^ and E"*" are functions of the three mutually commuting dimensionless operators: 7 °, 
zsj_ 7 ^ 7 ^, and (tt - 7 )^/^, where s_l = sign(ei?) and I = ^yhc/\eB±^\ is the magnetic length. In principle, they could also 
depend on the quasiparticle energy w. Taking into account that ( 7 °)^ = I and (zs_l 7 ^ 7 ^)^ = the Dirac structure of 
the operator-valued functions F~^ and E’*' can be written in the following form: 


= f + l°g + + *S-L7°7^7^/, 

E+ = A -I- 7 °/i -I- zs_L7^7^/i -I- *s_l 7 ° 7 ^ 7 ^A, 


( 6 ) 

(7) 


where /, /, g, g, A, A, g,, and g are functions of only one operator, (tt • To a large degree the physical 

meaning of the corresponding operators should be clear from their Dirac structure and symmetry properties. In 
particular, the eigenvalues of the first four of them (/, /, g, and g) will play the role of generalized wave-function 
renormalizations, while the others will be playing the role of the generalized the MC and QBE order parameters, i.e., 
the Dirac (parity-even) and Haldane (parity-odd) masses (i.e., A and A) and chemical potentials (i.e., g and g). As 
we will see, such an interpretation is also supported by the role of the corresponding parameters in the dispersion 
relations of quasiparticles, see Eqs. da and d below. 

The three mutually commuting operators, zs± 7 ^ 7 ^, and (tt • 7 )^?^, allow a common basis of eigenstates, 
|soj S 12 , ?T-)- The corresponding eigenvalues are Sq = ±1, S 12 = ±1, and —2n = —{2N -|- I -I- s_lSi 2 ), respectively, 
where A = 0, 1 , 2 , ... is the quantum number of the orbital angular momentum and S 12 is the sign of the pseudospin 
projection. By making use of the complete set of eigenstates |so,si 2 ,n), one can derive a very convenient generalized 
Landau level representation (GLLR) of the (inverse) quasiparticle propagator ([5]). (Eor details of the derivation, see 
Appendix A in Ref. The final form of the inverse full propagator is given by 


G-^(a;;r,r') = e*‘*'('’’‘'')G-^(a;; r - r'). 


( 8 ) 


where $(r, r') = — sj_ ^ ) is the well-known Schwinger phase in an external magnetic field in the Landau gauge 

A = (0,Bx). The translation invariant part of the inverse GLLR propagator is given by 


-G 2 


G r) — —i “I S 0 Wl/r!,(C) + ('SOMn.cr + A„_cr) -I-5^°^I/„_i(^)] 


n —0 ( 7=±1 


T* ^2 (7 ■ '")(/Ti,cr SQgn,a)L^_i{^)\Fso,s 


(9) 


where ^ /{2P) and L“(^) are the Laguerre polynomials. (Here, by definition, (C) = L„(^) and L“^g(^) = 0.) In 

Eq. dH), we also used the following set of eigenstate projectors in the Dirac space: Pso,si 2 = 7 ( 1 +so 7°)(1 + si 2 A_l 7 ^ 7 ^), 
as well as the following shorthand notations for the linear combinations of the order parameters: 

-j- fjAyj, fn,a — fn F ^fnj 9n^a — gn F ^gnt (^9) 






An,a = A, 


where tr = S 0 S 12 . Note that the parameters /„, fn, gn, gn, An, An, gn, and gn are associated with the nth Landau 
level. They are obtained by calculating the eigenvalues of the corresponding operators, introduced on the right-hand 
side of Eqs. and ©• 
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Similarly, the full GLLR propagator takes the form 

G(a;;r,r') 


( 11 ) 


where the Schwinger phase is exactly the same as in Eq. ([5]) and the translation invariant part of the propagator reads 




n—0 so=±l (7=±1 


.vt 


+*7^(7-Oe- 

(UJ 


fn,(T SQgn,<T I „ 

_i_ ,, \2 TP2 (' so,so<y i 

+ Mn.o-j ^ri,(T J 


( 12 ) 


The explicit form of the Landau level energies En^a are given by 


Eo,cr — trAo.o- — Aq + crAo, (13) 

En,a = \l 2n (/2 - gl ,^) vl/P + A2 n > 1. (14) 

The corresponding quasiparticles energies are determined by the location of the poles in the propagator (TT^ . namely 
^o,cr — ^0,(7 4“ Eq^(j and — f^n,(T i Eji u for n ^ 1. 

Note that the expressions for the Landau level energies in Eqs. 0 and (fTdl) shed additional light on the physical 
meaning of the various dynamical parameters, used in the ansatz of the full quasiparticle propagator. In particu¬ 
lar, as we can see, the combination of the wave-function renormalization parameters fn,a ~ 9n a determines the 
renormalization of the quasiparticle velocity. Also, we see that the absolute value of A„_cr plays the role of a mass. 


C. Schwinger-Dyson (gap) equation with static screening 


In order to derive the GLLR form of the Schwinger-Dyson (gap) equation, we will start from the standard coordinate 
form of the gap equation in the random-phase approximation (RPA) and assume that the interaction is provided by 
the l ong -range Goulomb force subject to static screening effects. Such a consideration will amend the analysis of 
Ref. |48j| . where the corresponding gap equation was analyzed in the approximation without screening. Our goal here 
is to illuminate the qualitative and quantitative role played by screening. 

It is important to emphasize that we will use a mean-field method that ignores the fluctuations of order parameters. 
While this is a common approximation used in numerous theoretical studies, there is no rigorous justification of its 
validity. Eormally, the corresponding fluctuations can destroy any long-range order in 2D when T ^ 0 and, thus, 
prevent any spontaneous breakdown of continuous symmetries. We argue semirigorously that this may not be the case 
for finite-size, high-quality graphene devices fabricated on substrates. The substrate could make certain aspects of 
dynamics in graphene effectively three-dimensional and, thus, tame the dangerous fluctuations (at least on the length 
scales of typical devices) that would destroy the long-range order in an ideal 2D graphene. Admittedly, however, this 
issue requires a more careful investigation in the future studies. 

After taking into account the condition of overall charge neutrality in graphene, we arrive at the following gap 
equation for the translation invariant part of the quasiparticle propagator [48l |: 

G-i(w;r) = S-\uj-,r) + (15) 

where D(u;; r) is the time-like component of the gauge (photon) field propagator that contains all the information 
about the interaction. As stated before, we will use an approximation with Goulomb interaction that includes the 
effects of static screening. In essence, this is the instantaneous approximation Q that neglects the retardation of 
the interaction. One might argue that this is a reasonable approximation because charge carriers are much slower 
than the speed of light. (See, however. Refs, [sOl suggesting that dynamical screening could affect the dynamics 
quantitatively.) Therefore, we use an energy independent polarization function, i.e., n(a;,k) ~ 11(0, k), in order to 
model the effects of screening in the photon propagator. In momentum space, the latter reads 


D{uj,k) ft! D(0,k) 


i 

eo [k + n( 0 ,k)] ’ 


(16) 


where eg is the dielectric constant, associated with the substrate. When the value of eg is large, the underlying 
dynamics becomes weakly coupled and the approximation for the gap equation (IlSp should become reliable. This is 
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not always the case in graphene. However, we expect that this may be applicable in the case of the highest quality 
graphene devices, fabricated on a thin hexagonal boron nitride substrate [23|. 

A convenient explicit form of the polarization function in the presence of an external magnetic field was calculated 
in Refs. [d^. In the approximation that neglects the wave-function renormalization and the Dirac masses, it reads 


n(o,k) = 


■EE 


Qnuiy) 


cosh^f^ 47r/^ 

n=0 A=+ \ 2T ) 


E E 


AA'/ .nF{\Mn) - nF{\'Mn') 


n,n'—0 A,A' —± 
Xn^X'n' 


AM„ - A'M„ 


where M„ = ^2nv'p/l^, T is the temperature, and the explicit form of function Q^n'iu) given by [4; 

+ AA'd - 


{n< 1 )! |n-n>| 

(n> - 


(17) 


(18) 


Here, by definition, y = k^P/2, n> =max(n,n') and n< =min(n,n'). We note that neglecting the wave-function 
renormalization and masses in the calculation of the polarization function should provide a reasonable approximation 
to leading order in weak coupling. Moreover, this may work also at moderate coupling because the bulk of the 
polarization effects appear to be determined by the total number of filled Landau levels and not as much by the 
details of the quasiparticle dispersion relations. 

For the derivation of the GLLR form of the gap equations, we refer the reader to Appendix B in Ref. [d^. The 
final set of equations reads 


for n > 0, and 


OO 

/^n,cr /^e ^ ^ ^ (E'n'',cr l^n' (E'n',cr H“ l^n',cr^ 


n'—O 

,(7 


[Tlp(^fin',a E^n',(T) ,(T “t“ /^n',^)] m 

^ /T 


(19) 


OO 

l^n,(7 Me “1“ ^^n,(7 — ^ ^ ^n' — — '^F(E^n',(T H“ 


n' — l 


jy, ^(T E'n^o') ^F(E'n\CT H“ fJ^n',<7')] I* 7 

■^n',(T ' 




00 ( 1 ) r 

\ ' ^n' —l,n—1 In',a 


2 , _ ^ Eji' (j 

n' — l ' 


F ^(T E'n^jCr) ^F(E'n^,cr “1“ l^n' ; 


9rL,a — 


aci 


00 ^( 1 ) 

E " [nF[yn',a - En',cr) “ nF{En',a + fJ-n',a)] , 


2 n E, 

n' — l 


( 20 ) 

( 21 ) 

( 22 ) 


for n > 1. In these equations, we introduced a dimensionless coupling constant, a = je^VF ~ 2.2, which is an 
analog of the fine structure constant for suspended graphene. 

The effect of screening in the above set of gap equations is implicit. It comes only through the modified values of 
the kernel coefficients Ida. 


K 


(p) 

m,n 


7” dk klF^]n{kl) 

Jq 2 f k + n( 0 , k) ’ 


with p = 0,1. 


(23) 


It is straightforward to calculate numerically the static polarization function n(0, k) as a function of dimensionless 
variable kl. We checked that, for sufficiently small and sufficiently large values of kl, the numerical results approach 
the expected asymptotical behavior at fc —?> 0 and k —>■ oo, respectively, derived in Eqs. (19) and (22) in Ref. (d^ . 
Then, by making use of the definition in Eq. (1231) . we obtain the numerical values of the kernel coefficients Km]n (with 
p = 0,1). The results for the first few Landau levels are presented in Tables ID and HIl respectively. By comparing these 
with the kernel coefficients in Ref. [d^, we observe that the screening effects substantially (i.e., by about a factor of 
2 ) decrease the values of Km]n- As we will see in the next section, this change strongly affects the magnitude of the 
order parameters (i.e., the symmetry breaking Dirac masses and chemical potentials) in the QH states with fractional 
filling of Landau levels. 
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TABLE 1. The kernel coefficients Kmln with the effects of static screening included. 


^m,Ti 

m = 0 

m = 1 

m = 2 

m = 3 

m = 4 

m = 5 

m = 6 

m = 7 

m = 8 

m = 9 

m = 10 

0 

0.0744 

0.0235 

0.0164 

0.0136 

0.0120 

0.0109 

0.0101 

0.0094 

0.0088 

0.0083 

0.0079 

1 

0.0235 

0.0602 

0.0222 

0.0156 

0.0129 

0.0114 

0.0104 

0.0096 

0.0090 

0.0085 

0.0081 

2 

0.0164 

0.0222 

0.0531 

0.0211 

0.0149 

0.0124 

0.0109 

0.0100 

0.0093 

0.0087 

0.0083 

3 

0.0136 

0.0156 

0.0211 

0.0486 

0.0202 

0.0144 

0.0120 

0.0106 

0.0097 

0.0090 

0.0085 

4 

0.0120 

0.0129 

0.0149 

0.0202 

0.0453 

0.0195 

0.0140 

0.0116 

0.0103 

0.0094 

0.0088 

5 

0.0109 

0.0114 

0.0124 

0.0144 

0.0195 

0.0428 

0.0189 

0.0137 

0.0113 

0.0100 

0.0092 

6 

0.0101 

0.0104 

0.0109 

0.0120 

0.0140 

0.0189 

0.0407 

0.0184 

0.0133 

0.0111 

0.0098 

7 

0.0094 

0.0096 

0.0100 

0.0106 

0.0116 

0.0137 

0.0184 

0.0390 

0.0180 

0.0131 

0.0109 

8 

0.0088 

0.0090 

0.0093 

0.0097 

0.0103 

0.0113 

0.0133 

0.0180 

0.0375 

0.0176 

0.0128 

9 

0.0083 

0.0085 

0.0087 

0.0090 

0.0094 

0.0100 

0.0111 

0.0131 

0.0176 

0.0362 

0.0172 

10 

0.0079 

0.0081 

0.0083 

0.0085 

0.0088 

0.0092 

0.0098 

0.0109 

0.0128 

0.0172 

0.0351 


TABLE 11. The kernel coefficients Km]ri with the effects of static screening included. 



m = 0 

m = 1 

m = 2 

m = 3 

m = 4 

m = 5 

m = 6 

m = 7 

m = 8 

m = 9 

m = 10 

0 

0.0509 

0.0142 

0.0084 

0.0064 

0.0055 

0.0051 

0.0048 

0.0046 

0.0043 

0.0041 

0.0039 

1 

0.0142 

0.0901 

0.0283 

0.0172 

0.0131 

0.0111 

0.0099 

0.0093 

0.0088 

0.0084 

0.0081 

2 

0.0084 

0.0283 

0.1244 

0.0419 

0.0260 

0.0198 

0.0167 

0.0148 

0.0137 

0.0129 

0.0123 

3 

0.0064 

0.0172 

0.0419 

0.1553 

0.0550 

0.0346 

0.0264 

0.0222 

0.0197 

0.0181 

0.0169 

4 

0.0055 

0.0131 

0.0260 

0.0550 

0.1839 

0.0677 

0.0431 

0.0330 

0.0276 

0.0244 

0.0224 

5 

0.0051 

0.0111 

0.0198 

0.0346 

0.0677 

0.2107 

0.0801 

0.0515 

0.0395 

0.0331 

0.0292 

6 

0.0048 

0.0099 

0.0167 

0.0264 

0.0431 

0.0801 

0.2361 

0.0922 

0.0597 

0.0459 

0.0385 

7 

0.0046 

0.0093 

0.0148 

0.0222 

0.0330 

0.0515 

0.0922 

0.2602 

0.1039 

0.0678 

0.0523 

8 

0.0043 

0.0088 

0.0137 

0.0197 

0.0276 

0.0395 

0.0597 

0.1039 

0.2834 

0.1154 

0.0758 

9 

0.0041 

0.0084 

0.0129 

0.0181 

0.0244 

0.0331 

0.0459 

0.0678 

0.1154 

0.3056 

0.1266 

10 

0.0039 

0.0081 

0.0123 

0.0169 

0.0224 

0.0292 

0.0385 

0.0523 

0.0758 

0.1266 

0.3271 


III. NUMERICAL RESULTS 


In this section, we present our numerical solutions to the GLLR gap equations (1191) through (0^ in the model 
with static screening effects. In the calculation, we use the Newtonian iteration algorithm and replace the original 
infinite set of gap equations with a truncated set of Umax = 50 equations. In accordance with such a truncation, we 
also impose a sharp cutoff in the summation over the Landau level index at rimax = 50. (The numerical tests reveal 
that the cutoff at rimax = 100 does not affect any qualitative features of the solutions and only slightly changes the 
numerical results for the dynamical parameters in the first few Landau levels.) 

Taking into account that each Landau level (n < rimax) has two possible spin states and each of them is characterized 
by eight different dynamical parameters, see Eq. OT, we have a total of IGrimax independent parameters that should 
be determined by solving lOrimax coupled algebraic equations, see Eq. (fTOl) through (1^ . [Strictly speaking the number 
of independent parameters is Ibrimax + 4, where the additional four corresponds to the n = 0 Landau level, which is 
special.] 

Before proceeding with the numerical analysis, it is useful to note that the coupled sets of gap equations for the 
two spin states, s =4,,t) have exactly the same form. They differ only in the value of the chemical potential, i.e., 
fjL^ = — Z and /rj, = + Z, where Z = ^bB is the Zeeman energy. We also note that the gap equations are the 

same for cr = ±1 when the dynamical parameters fn,a, gn,cT, gn,a and crAn^a are treated as independent variables. 
These two observations allow us to greatly reduce the numerical cost of calculations. We first consider the generic 
set of equations in a finite range of chemical potentials and tabulate possible solutions for the above-mentioned four 
independent parameters. Usually, there exist multiple solutions that differ in their symmetry properties. The final 
solutions for the complete set of dynamical parameters with fixed values of s =4,,t cr = ±1 are obtained by 
combining the distinct generic solutions for properly shifted values of the chemical potential. 












8 


TABLE III. Values of the wave-function renormalization /„ for several values of the chemical potentials with screening effects 
considered, compared with the solutions in Ref. with screening effects neglected in the parentheses. (Note that the 
corresponding data in Ref. had typos that were corrected here.) 




/2 

h 

U 

h 

h 

Pe < V2ei 

1.084 (1.270) 

1.072 (1.243) 

1.065 (1.227) 

1.060 (1.214) 

1.057 (1.205) 

1.054 (1.197) 

\/2ti < |pe| < Vde; 

1.045 (1.194) 

1.066 (1.224) 

1.063 (1.217) 

1.059 (1.208) 

1.056 (1.201) 

1.054 (1.194) 

y/^ei < pe < \/Qti 

1.037 (1.166) 

1.042 (1.177) 

1.058 (1.200) 

1.057 (1.199) 

1.055 (1.194) 

1.052 (1.189) 

%/6ei < pe < -s/Sq 

1.033 (1.150) 

1.035 (1.156) 

1.039 (1.165) 

1.052 (1.184) 

1.052 (1.185) 

1.051 (1.182) 

%/8ei < |/ie| < VWei 

1.031 (1.138) 

1.032 (1.142) 

1.034 (1.148) 

1.037 (1.156) 

1.048 (1.172) 

1.049 (1.174) 

VWei < pe < •s/l2ei 

1.029 (1.128) 

1.030 (1.132) 

1.031 (1.136) 

1.032 (1.141) 

1.035 (1.148) 

1.045 (1.162) 

%/l2ei < |pe| < Vide; 

1.027 (1.121) 

1.028 (1.123) 

1.028 (1.127) 

1.029 (1.130) 

1.031 (1.135) 

1.034 (1.141) 


In the numerical analysis, it is convenient to express all physical quantities in units of the Landau energy scale, 

ei = ^Jhvl,\eBjL\/c « [T] meV Ri SOOa/B^ [T] K, (24) 

where the value of the magnetic field is measured in teslas. It may be appropriate to emphasize here that, while e/ 
is determined by the component of the magnetic field orthogonal to the plane of graphene Bj^, the Zeeman energy 
Z = ^bB ~ 5.8 X [T] meV ~ 0.67i? [T] K is proportional to the absolute value of B. Despite the fact that 

B± < B, the Zeeman energy is generically much smaller than the Landau energy scale (|24)) . This changes only 
in the case when the magnetic field becomes nearly parallel to the plane of graphene, i.e., when B± <C B. In the 
analysis below, we do not consider such a limiting case. In general, the dimensionless Zeeman energy is given by 
z = Z/ei = 2.2 X 10~^y/B [T]/ cosOb, where 6b is the angle of the magnetic field tilt. In the numerical analysis 
below, we will fix the value of the dimensionless Zeeman energy to be z = 0.015. This formally corresponds to 
5/cos 03 a: 46.5 T. 


A. Fermi velocity renormalization in weak magnetic field 

In a weak magnetic field, the effect of symmetry breaking dynamical parameters is expected to be negligible. 
However, even in this regime it is interesting to explore the implications of the long-range Coulomb interaction on the 
quantum Hall effect in graphene. In particular, the interaction is expected to renormalize the Fermi velocity, which can 
be extracted from the dynamically modified expressions for the Landau level energies [i^ . In absence of QHF and MC 
order parameters, as we see from Eq. o, the energies are given by En = fnVpV^^^i, implying that the renormalized 
Fermi velocity is determined by the wave-function renormalization parameters /„, namely = fnVp- Without 

screening of the Coulomb interaction, the numerical values for the wave-function renormalization parameters /„ were 
previously reported in Ref. [i^ . The corresponding results with the effects of screening were reported in Ref. [s^ for 
the QH states with filling factors p = + 1/2). In the latter study, in fact, the dynamical parameters and A„ 

were also properly accounted for. By noting that the values of /„ depend not only on the chemical potential, but also 
on the Landau level index n, we conclude that the same is also true for the renormalized Fermi velocity. This is a 
very interesting theoretical prediction that could be easily tested in o ptic al experiments, for example, via a systematic 
study of absorption/transmission lines for each of the QH states |5ll - l^ . 

In the context of the optical transitions, the effect of interactions can be conveniently quantified by measuring the 
deviations of the measured energies of transitions from the free theory predictions [52| , 

^En,n' = En' ± En = ± v^)e/ -|- aeiCn,n', (25) 

where the minus sign corresponds to transitions between states with negative energies and states with positive energies. 
Strictly speaking, the above definition of AEn^n> assumes transitions between states with the same spin. In the case 
of transitions between different spin states, an extra ±2Z correction should be added on the right-hand side. Here, 
the information about the wave-function renormalization is captured by the following set of dimensionless parameters 
Cn,n't i.e.. 


Cn,n' = - I) ± ^(/n - 1). 


a 


(26) 
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TABLE IV. Parameters in the QH state with the filling factor u = 2 m the model with static screening. The values in 

parentheses are obtained in an approximation with the wave-function renormalization effects included, but all QHF and MC 
order parameters neglected. 


^n,n' 

n' = 1 

n' = 2 

n' = 3 

n' = 4 

n = 0 

0.082 (0.054) 

0.093 (0.065) 

0.099 (0.072) 

0.103 (0.077) 

n — —1 

0.100 (0.108) 

0.111 (0.119) 

0.117 (0.126) 

0.122 (0.131) 

n = —2 

0.108 (0.119) 

0.119 (0.130) 

0.125 (0.138) 

0.129 (0.143) 

n = —3 

0.113 (0.126) 

0.124 (0.138) 

0.130 (0.144) 

0.134 (0.150) 

n = —4 

0.116 (0.131) 

0.127 (0.143) 

0.133 (0.150) 

0.138 (0.220) 


As we see, in absence of symmetry breaking, which is the case in weak magnetic fields, the values of parameters Cn,n' 
are directly related to the quasiparticle velocity renormalizations. As we claim here, for each QH state, the corre¬ 
sponding renormalizations are functions of the Landau level index n. In experiment, the complete set of parameters 
/„ with n > 1 could be extracted by measuring the transition energies between the lowest Landau level (which is free 
from from the corresponding renormalization effect) and higher Landau levels. The values of /„, extracted in this 
way, would be sufficient to calculate the values of Cn,n' for transitions between various higher Landau levels. The 
latter could be also compared to the actual measurements and, in the case of agreement, one would have a nontrivial 
test of the self-consistency of the GLLR ansatz used here. From a theoretical point of view, it will be perhaps even 
more interesting if deviations from the relations in Eq. (1261) are observed. 

In a general case with nonvanishing QHF and MC order parameters, the expressions for Cn,n' parameters should be 
corrected because the Landau level energies are modified, see Eqs. (US and (HI- The magnitude of the corresponding 
corrections is of the order of \^in,a — Mel and for Landau levels n > 1. (Note that the correction due to 

the QHF order parameter should be interpret as part of the energy measured with respect to the thermodynamical 
potential /ig.) In the case of transitions to/from the lowest Landau level (n = 0), the corrections are ±crAo,cr, see 
Eq. (1131) . Below, we present our numerical results for the parameters in several QH states in both approximations, 
i.e., with and without inclusion of the QHF and MC order parameters. 

By utilizing the same value of the cutoff, Umax = 100, in the numerical calculation as in Ref. (48j| . but also 
including the effects of screening, we straightforwardly obtain the wave-function renormalization parameters /„. The 
corresponding zero temperature results for the first few Landau levels are presented in Table IHII We note that they 
differ slightly from those in Ref. [s^ because here all dynamical parameters such as /i„ and A„ were neglected. For 
comparison, in Table m we also list in the parentheses the numerical results from Ref. [i^, obtained without the 
screening effects. Different rows in Table IHII correspond to different choices of the chemical potential in the gaps 
between a completely filled nth Landau level and a completely empty (n -I- l)th Landau level. 

By comparing the results with and without screening in Table Hill we find that the effect of wave function renormal¬ 
ization goe s from about 14% to 27% down to about 3% to 8%, which appears to be even smaller than the prediction 
in Ref. |54| . It is curious to explore in detail if the logarithmic running of the coupling constant could explain such 
a difference. By making use of the definition in Eq. (1261) and our numerical results for various QH states, we readily 
calculate Cn,n' parameters. 

The representative results for Cn,n' parameters are presented in Tables 113 and El for the QH states with the 
filling factors v = 2 and v = 6. Because of the particle-hole symmetry, the results for the corresponding negative 
filling factors can be obtained as follows: Cn,n'{—^) = C-r»' These results appear to be about two to four 

times smaller than the results in the absence of screening [^. In calculation, we took into account the effect of 
nonvanishing QHF and MC order parameters. For comparison, in parenthesis we also show the results for the same 
parameters in the approximation with the QHF and MC order parameters neglected. It appears that the role of such 
parameters is not negligible. 


B. Quantum Hall states in strong magnetic field 

In this section, we study QH states with different integer filling factors at zero temperature. We will start by first 
considering the states associated with the n = 0 and n = 1 Landau levels. We will show, in particular, that the 
qualitative features of the phase diagram obtained in Ref. remain qualitatively the same after the inclusion of the 
static screening. At the same time, the quantitative changes will be substantial. 

As stated earlier, the gap equations (fTOl) through (1^ allow a large number of solutions with different types of 
Haldane/Dirac masses (A„, A„) and chemical potentials (/r„, (i„). In order to identify the true ground state among 
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TABLE V. Parameters Cn^n' in the QH state with the hlling factor i/ = 6 in the model with static screening. The values in 
parentheses are obtained in an approximation with the wave-function renormalization effects included, but all QHF and MC 
order parameters neglected. 


^n,n' 

n' = 2 

n' = 3 

n' = 4 

n = 1 

0.054 (0.031) 

0.066 (0.041) 

0.073 (0.047) 

n = 0 

0.088 (0.060) 

0.100 (0.070) 

0.107 (0.076) 

n = —1 

0.103 (0.089) 

0.115 (0.100) 

0.122 (0.105) 

n = —2 

0.109 (0.121) 

0.121 (0.131) 

0.128 (0.136) 

n = —3 

0.113 (0.131) 

0.125 (0.141) 

0.132 (0.147) 

n = —4 

0.116 (0.136) 

0.128 (0.147) 

0.135 (0.152) 


them, we compare their free energies. (For the explicit expression of the free energy, see Appendix C in Ref. (d^.i 
In the model at hand, the choice of the corresponding states is strongly affected by the Zeeman interaction energy, 
which is one of the main factors in driving the vacuum alignment. Taking into account that there may exist a large 
number of other symmetry breaking effects, e.g., various on-site repulsion interaction terms [l^ . the actual nature of 
the ground states should be accepted with caution. Nevertheless, the study below is quite informative: it reveals the 
quantitative effects of the static screening and role of the running of the dynamical parameters in the QH regime of 
graphene. 

From the symmetry viewpoint, there are two types of Dirac masses and chemical potentials for quasiparticles of 
each spin orientation, s =4,,t- The parameters of the first type (i.e., A„ and /r„) are singlets, while the parameters of 
the second type (i.e., A„ and jin) are triplets with respect to the flavor Us(2) subgroups. It is natural to expect that 
the states with different symmetry properties compete. However, it should be emphasized that the corresponding 
competition is not necessarily between the MC and QHF scenarios because both types of order parameters may belong 
to the same representations of the flavor symmetry. In fact, as our results show, the two types of order parameters 
generically coexist in all QH states. (This was also emphasized in Ref. (d^.l 

From general considerations based on the structure of the Landau levels, it is expected that there exist at least four 
different classes of QH states with the following series of filling factors: (i) i/ = An + 2, (ii) v = An (with a possible 
exclusion of the rather unique iz = 0 state in a class of its own), (iii) v = An — 1, and (iv) v = An + 1. 

The first series of states with the filling factors v = An + 2 describes the “normal” QH states with the complete 
filling of the (nearly) degenerate Landau levels. They do not have or require symmetry breaking and are resolved even 
at relatively weak magnetic fields. The simplest realization of the v = An states could be provided by a dynamically 
enhanced Zeeman splitting of Landau levels. In this case, the ground states have Ui{2) x U^{2) symmetry. While 
this is also the prediction of the model at hand, we should emphasize that other realizations of the v = An QH states 
are possible. In fact, the n = 0 state, which formally belongs to this series, is likely to have a different origin [dj. 
The remaining two series of states with the filling factors n = An Az I are less controversial. They require spontaneous 
symmetry breaking at least down to t/ 4 .(l) x U^{2) or Ui{2) x 

In our analysis at sufficiently small values of the chemical potential, we find a number of different solutions, 
associated with the lowest Landau level (n = 0) and integer filling factors n = 0, ±1, ±2. The order of appearance and 
the competition of different types of solutions appear to be the same as in Ref. [^, see Figure 3 there. By taking into 
account the screening of the Coulomb interaction, however, we find that the actual values of dynamical parameters 
change substantially. The corresponding results are listed in Table IVII For comparison, in parenthesis we also list 
the previous results in the model without screening [d^. As we see, the effect of screening is to suppress the relevant 
dynamical parameters by about a factor of 3. The same is true for the magnitude of the energy gaps in the states 
with the filling factors z/ = 0, ±1. 

The analysis of the QH states, associated with the n = 1 Landau level, is done in the same way. Here again, the 
order of appearance and the competition of different types of solutions appear to be exactly the same as in Ref. (d^ . 
see Figure d there. The actual values of the dynamical parameters in the model with screening are listed in Table IVTIl 
The corresponding results in the model without screening are given in parenthesis. By comparing the two sets of 
data, we find that screening leads to a suppression of the relevant dynamical parameters by a factor of 3 to 5. 

Because of the long-range nature of the Coulomb interaction, the dynamical QHF (/x„, jin) and MC (A„, A„) 
order parameters are nontrivial functions of the Landau level index n. The corresponding “running” of the dynamical 
parameters is an important feature of the model at hand. From theoretical viewpoint, it is essential to provide a 
realistic description of the low-energy dynamics in graphene, where the role of order parameters diminishes with 
increasing the quasiparticle energy. This is in contrast to a common mean-field analysis of models with pointlike 
interactions, where the order parameters affect either (i) only the nearest filled Landau level or (ii) all levels in the 
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TABLE VI. Gap parameters for the solutions when the Fermi energy close to the lowest Landau level. The solutions for 
JT-max = 50 with Screening effects considered are compared with the solutions in Ref. for rimax = 5 neglecting screening 
effects in the parentheses. 


V 


X eff 

^o,t 



X efF 



Ao,t 



A04 

-2 


0.000 (0.000) 


0.000 (0.000) 


0.082 (0.227) 


0.082 (0.227) 

-1 


0.000 (0.000) 


0.082 (0.227) 


0.082 (0.227) 


0.000 (0.000) 

0 


0.000 (0.000) 


0.000 (0.000) 


0.082 (0.227) 


-0.082 (-0.227) 

1 


0.082 (0.227) 


0.000 (0.000) 


0.000 (0.000) 


-0.082 (-0.227) 

2 


0.000 (0.000) 


0.000 (0.000) 


-0.082 (-0.227) 


-0.082 (-0.227) 

TABLE VII. Gap parameters for the solutions when the Fermi energy close to the Landau level at n 

= 1. The solutions for 

^max — 

50 with 

screening effects considered are 

compared with the solutions in 

Ref. [lil for 

n-max — 

5 neglecting screening 

effects in the parentheses. 











Spin 

X efF 

^0 

Ao 

/i 

pi — Ps 

Ai 

Ai 

/2 

P 2 — Ps 

A2 

A2 


t 

0.000 

-0.082 

1.078 

0.013 

-0.014 

0.000 

1.065 

0.009 

-0.010 

0.000 


(0.000) 

(-0.227) 

(1.143) 

(0.053) 

(-0.068) 

(0.000) 

(1.112) 

(0.040) 

(-0.052) 

0.000 


4 

-0.013 

-0.095 

1.058 

0.050 

-0.010 

0.004 

1.062 

0.022 

-0.009 

0.000 


(-0.051) 

(-0.278) 

(1.105) 

(0.148) 

(-0.049) 

(0.018) 

(1.102) 

(0.091) 

(-0.046) 

(0.006) 

u — 4 

t 

0.000 

-0.082 

1.078 

0.013 

-0.014 

0.000 

1.065 

0.009 

-0.010 

0.000 


(0.000) 

( -0.227) 

(1.143) 

(0.053) 

(-0.068) 

(0.000) 

(1.112) 

(0.040) 

(-0.052) 

(0.000) 


i 

0.000 

-0.108 

1.038 

0.087 

-0.006 

0.000 

1.059 

0.034 

-0.009 

0.000 


(0.000) 

(-0.330) 

(1.066) 

(0.243) 

(-0.031) 

(0.000) 

(1.093) 

(0.143) 

(-0.039) 

(0.000) 

p — 5 

t 

-0.013 

-0.095 

1.058 

0.049 

-0.010 

0.004 

1.062 

0.021 

-0.009 

0.000 


(-0.051) 

(-0.278) 

(1.105) 

(0.148) 

(-0.049) 

(0.018) 

(1.102) 

(0.091) 

(-0.046) 

(0.006) 


1 

0.000 

-0.108 

1.038 

0.087 

-0.006 

0.000 

1.059 

0.034 

-0.009 

0.000 


4- 

(0.000) 

(-0.278) 

(1.066) 

(0.245) 

(-0.031) 

(0.000) 

(1.093) 

(0.143) 

(-0.039) 

(0.000) 

u — 6 

t 

0.000 

-0.108 

1.038 

0.087 

-0.006 

0.000 

1.059 

0.034 

-0.009 

0.000 


(0.000) 

(-0.330) 

(1.066) 

(0.243) 

(-0.031) 

(0.000) 

(1.093) 

(0.142) 

(-0.039) 

(0.000) 


4 

0.000 

-0.108 

1.038 

0.087 

-0.006 

0.000 

1.059 

0.034 

-0.009 

0.000 


(0.000) 

(-0.331) 

(1.066) 

(0.245) 

(-0.031) 

(0.000) 

(1.093) 

(0.143) 

(-0.039) 

(0.000) 


same way. 

Our numerical results for the running dynamical parameters fn,s, ^■n,s, ^n,s, An.s, and are shown in Fig. [1] 
through Fig. |4l In line with our discussion of the four different classes of QH states with different filling factors, we 
show the results for = 4n + 2, i' = 4n, i/ = 4n — 1, and = 4n + 1 states in separate panels. They are characterized 
by different ground state symmetries. In order to avoid overcrowding the hgures, we showed the results only for a few 
states in the series, that correspond to partially or fully filled n = 0, n = 1, and n = 2 Landau levels. It is natural to 
expect that the other states in the same series, associated with filling of higher Landau levels, share essentially the 
same qualitative features. 

In Figs, m - m the results for the spin-up and spin-down quasiparticle states are represented by the same types 
of filled and unfilled symbols, respectively. The universal property of all dynamical parameters is that their values 
approach the free model limit at large n. Additionally, we see that often the running parameters acquire their largest 
absolute values in the Landau levels near the Fermi energy. These features were expected, of course. 

The results for the wave-function renormalization parameters fn,s as functions of n are shown in Fig. [TJ As we see, 
the results are qualitatively the same for all four different classes of QH states. This suggests that the wave-function 
renormalization parameters and, thus, the renormalized Fermi velocities are largely determined by the long-range 
(screened) Coulomb interaction and not very sensitive to the effects of the dynamically generated symmetry breaking 
terms. 

In the states with the even filling factors v = 4n and zz = 4n -|- 2, only the singlet order parameters 
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FIG. 1. (Color online) The wave-function renormalization parameters fn,s as functions of the Landau level index n for four 
different types of QH states with filling factors iz = 4n — 1, v = 4n, v — An + 1, and i/ = An + 2. 


and An,s are generated. In both cases, the dynamical parameters have comparable magnitudes and the ground state 
symmetry is C/f (2) x 174,(2). While the role of nontrivial fin,i — ^n,s is critical in the series of states with 

the filling factors v = An, it is not the case for the states with iz = 4n -|- 2. Indeed, in the n = An states, it is the 
singlet order parameters that determine the energy gaps of the QH states. They are critical because the corresponding 
dynamical energy gaps are generated to be much larger than the bare Zeeman splitting. This is in contrast to the 
ly = An+ 2 states, which are characterized by rather large gaps (of order ei) between Landau levels, which cannot be 
affected substantially by relatively small corrections due to fin,i — ^n,s- 

Our results for the v = 4n, associated with higher Landau levels, appear to be in qualitative agreement with the 
experimental results reported in Refs. [2^, [2^. The spin-polarized nature of the corresponding states is supported by 
the observed increase of the activation gaps as functions of the in-plane component of the magnetic field. 

In agreement with our general symmetry arguments, the triplet order parameters fln,s a-nd A„ ,, are generated only 
in the states with the odd filling factors zz = 4n ± 1, see Fig. ID We should note, however, that in both types of states 
the singlet order parameters — fJ-n.A and A„_s are generated as well, see Figs. [5] and [3] This is not surprising 
since they do not break any additional symmetries. The qualitative difference between the states with zz = 4n — 1 
and zz = 4n -|- 1 is that the triplet order parameters are generated either exclusively for the spin-down quasiparticles 
(zz = 4n — 1) or exclusively for the spin-up quasiparticles (zz = 4n -4- 1). It is not immediately clear whether these 
results are in perfect agreement with the experimental data in Ref. [23|. Our theoretical model appears to capture 
some of the key qualitative features of the experimental data. For example, the corresponding states are characterized 
by a spontaneous breakdown of the flavor symmetry and the gaps are not particularly sensitive to the Zeeman energy. 
On the other hand, the exact symmetry breaking pattern in the observed states with filling factors zz = 4n ± 1 is not 
clear. In theory, the ground state symmetry is either Ui{l) x U^{2) or Ui{2) x t7|'(I). While this does not contradict 
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FIG. 2. (Color online) The chemical potential differences fin,a — fia as functions of the Landau level index n for four different 
types of QH states with filling factors v = 4n — 1, v — 4n, u = An + 1, and ;/ = 4n + 2. 


the measurements, the actual symmetry could in principle be lower, i.e., t/ 4 ,(l) x tl-t-(l). One way to resolve the issue 
unambiguously is to investigate the spectrum of quasiparticles in detail, e.g., by studying systematically all allowed 
optical transitions. We hope that this will be done in future investigations. 


C. Temperature dependence of the energy gaps 


As is clear from our discussion in the previous subsection, the QH states with the filling factors i/ = 4n and 
11 = An ± I are characterized by the energy gaps that are largely determined by the dynamically generated QHF and 
MC order parameters. (Recall that the energy gaps in the ^ = 4n + 2 states are of the order of the Landau energy 
scale ei and, thus, are not affected much by the dynamical order parameters.) The corresponding parameters and, 
therefore, the energy gaps are expected to have a strong temperature dependence. In this subsection, we study such 
a dependence in detail. 

The numerical analysis of the GLLR gap equations (fT^ through (l22l) at nonzero temperature is qualitatively the 
same as at T = 0. After determining the ground states at various filling factors as a function of temperature, we 
can straightforwardly extract the temperature dependence of the energy gaps. For different types of the QH states, 
associated with various fillings of the lowest three Landau levels (n = 0,1, 2), our numerical results are shown in Figs. [5] 
and [6] There are several universal features of these results: (i) the gaps decrease monotonically with temperature, 
(ii) the overall size of the gaps decreases with increasing the Landau level index n. We also find that, for a fixed n, 
the gap functions in the states with v = An — 1 and ii = An+ 1 are almost exactly the same, see Fig. [S] 

By comparing the results in Fig. [5] with those in Fig. [SJ we see that the temperature dependence of the energy 
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FIG. 3. (Color online) The singlet Haldane masses as functions of the Landau level index n for four different types of QH 
states with filling factors v — An — 1, i/ = An, v — An + 1, and u — An + 2. 


gaps in the QH states with filling factors z/ = 4n ± 1 is qualitatively different from that in the ly = An states. By 
taking into account the very different symmetry properties of the corresponding ground states, this is not surprising 
at all. The = 4n ± 1 states are characterized by the triplet order parameters jln,s and which vanish in a 

symmetry restoring phase transition at the critical temperature Tf., see Fig. [5] The transition appears to be either a 
second-order, or a weakly first-order transition. The temperature dependence of the energy gaps in the z/ = 4rz ± 1 
states is fitted quite well by the following function: 


AE = AE^°'> 1 



(27) 


where AE^^'> is the energy gap at zero temperature and Tc is the critical temperature. Note that the approximate 
numerical values of the zero temperature gaps are 0.165ez, O.lSOe; and 0.112ez for the i/ = An±l states, associated 
with the Landau levels n = 0,1,2, respectively. The corresponding approximate values of the critical temperatures 
are O.OdSe;, 0.033e/ and 0.029e;, respectively. 

In contrast, the v = An states are characterized by the singlet order parameters ^n,i — ^n,s, which have 

the same symmetry as the Zeeman term and, thus, remain nonzero even at large temperatures, see Fig.[S) As a result, 
the corresponding transitions from the low-temperature n = An states with a dynamically enhanced Zeeman splitting 
to the high-temperature states without such an enhancement are generically smooth crossovers. In the case of the 
Z 2 = 0 state, however, we find a sign of a small discontinuity in the temperature dependence of the energy gap around 
T ~ 0.04ez, suggesting a weak hrst-order transition. Of course, such a transition is not related to a restoration of any 
exact symmetry and, thus, can be viewed as accidental. 
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FIG. 4. (Color online) The triplet chemical potentials fin,a and the Dirac masses A„_s as functions of the Landau level index 
n for the QH states with the filling factors v = An ±1. 


IV. DISCUSSION 

In this paper, we utilized a highly flexible GLLR representation for the fermion propagator to describe the QH states 
with integer filling factors in the low-energy model of graphene with long-range Coulomb interaction. By including the 
static screening effects in an external magnetic field, we amended the earlier study of Ref. [i^. While using a similar 
mean-held approximation, we found that the static screening has a substantial suppression effect on the dynamical 
order parameters in all QH states with spontaneously broken symmetry. Also, the deviations of the wave-function 
renormalization from 1 and the dynamical corrections to the Fermi velocity became noticeably suppressed. 

By making use of the framework that naturally incorporates the running of the dynamical parameters with the 
Landau level index n, we observed that the largest absolute values of such parameters are typically obtained for the 
Landau level near the Fermi energy. In the limit of large n, all dynamical parameters approach the corresponding 
bare values. By making use of numerical methods, we studied in detail the behavior of the dynamical parameters in 
all four different types of the QH states with the filling factors i/ = 4n -I- 2, v = 4n, and v = AnAz \ that have different 
symmetry properties. At weak fields, we also analyze the running of the renormalized Fermi velocity with the Landau 
level index n. 

In the low-energy model used, the states with the filling factors v = An + 2 and v = An have the C/t(2) x 174.(2) 
symmetry. They are characterized by the singlet type of the QHF and MC order parameters with respect to C/s (2) for 
both s =4- and s =\. While the role of the corresponding singlet parameters in the v = An+ 2 states is negligible, it 
is profound in the v = An states, where they determine the magnitude of the dynamically enhanced Zeeman splitting. 
This finding qualitatively agrees with the data reported in Refs. [2^, [2^, which supports spontaneous symmetry 
breaking and spin polarization in v = An states (with the exception of the zz = 0 state). The states with the odd 
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V^i T/ei 

FIG. 5. (Color online) The energy gaps as functions of temperature for the QH states with the filling factors ty = 4n dz 1. 



T/c, 


FIG. 6. (Color online) The energy gaps as functions of temperature for the QH states with the filling factors ly = 4n. 


filling factors ly = 4n iL 1 have rather different properties and are characterized by a lower ground state symmetry, 
i.e., either [/'|-(1) x 17^(2) or Uf(2) x f7j,(l). In such states, in addition to the singlet order parameters, there are two 
types of spontaneously generated triplet parameters. The latter play a critical role in the realization of the ly = 4niLl 
states by reducing the symmetry so as to allow the appropriate lifting of the fourfold degeneracy, as well as needed 
partial filling of Landau levels. These properties appear to be in agreement with the experimental data in Ref. [2^, 
even though the data was not sufficient to establish unambiguously the underlying symmetry breaking pattern. 

By extending the analysis to the case of nonzero temperature, we studied the energy gaps in the QH states with the 
filling factors ly = 4n and ^ = 4n ± 1. We found that the symmetry breaking (triplet) order parameters, describing 
the ly = 4n iL 1 states, vanish at a certain critical value of temperature Tc, where a second-order or weakly first- 
order transition occurs. In terms of the zero temperature gap the results for the critical temperature are 

approximately given by Tc Ri 0.26Ai?(°) for all jz = 4n ± 1 states associated with filling of the first few Landau 
levels. On the other hand, we do not detect a well-defined symmetry-restoring phase transition in the temperature 
dependence of the energy gaps in the iz = 4n states. Instead, we find a smooth crossover, which is consistent with the 
fact that the corresponding states have no real symmetry-breaking order parameters. 

The results of this study unambiguously suggest that the QH states of graphene with various integer filling factors 
are generically characterized by a rather large number of dynamical parameters that have a nontrivial running with 
the Landau level index n. The corresponding rich structure of Landau levels in the QH regime could be probed 
in detail via a systematic study of transition lines in optical experiments such as those reported in Refs. |5ll - l53l| . 
The temperature dependence of the energy gaps obtained in this study appears to be in a qualitative agreement 
with the measurements of the activation energies in Ref. m, El) EE EE, [20, EE Ell ■ while such an agreement is 
encouraging, it remains to be seen whether the theoretical predictions of the GLLR formalism could be matched 
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quantitatively to the experimental data. Before such a comparison is attempted, however, one may need to reevaluate 
the model approximations used in the present study. In particular, one should address the precise quantitative role of 
(i) dynamical screening, (ii) fluctuations of the order parameter fluctuations, and (iii) nonzero quasiparticles widths. 
Hopefully, all such effects in GLLR formalism will be addressed in the future studies. 
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